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ABSTRACT 

We report on the serendipitous discovery of a new X-ray source, SAX J1802.7- 
2017, ~ 22' away from the bright X-ray source GX 9+1, during a BeppoSAX 
observation of the latter source on 2001 September 16-20. SAX J1802. 7-2017 
remained undetected in the first 50 ks of observation; the source count rate in 
the following ~ 300 ks ranged between 0.04 c/s and 0.28 c/s, corresponding to 
an averaged 0.1-10 keV flux of 3.6 x 10~^^ ergs cm~^ s~^. We performed a timing 
analysis and found that SAX J1802. 7-2017 has a pulse period of 139.612 s, a 
projected semimajor axis of a^smi ~ 70 It-s, an orbital period of ~ 4.6 days 
and a mass function f{M) ~ 17 ± 5 Mq. The new source is thus an accreting 
X-ray pulsar in a (possibly eclipsing) high mass X-ray binary. The source was 
not detected by previous X-ray astronomy satellites, indicating that it is likely a 
transient system. 

Subject headings: stars: neutron — stars: magnetic fields — pulsars: general — 
pulsars: individual: SAX J1802. 7-2017 — X-ray: binaries 



1. Introduction 

High mass X-ray binaries (hereafter HMXBs) are young systems, and the neutron stars 
(NSs) that are often hosted in them usually have a strong magnetic field {B ~ 10^^ Gauss). 
Accretion onto these NSs occurs via capture of stellar wind matter and Roche-lobe overflow. 
HMXBs comprise two main subgroups: a) the supergiant systems, and b) the Be star systems. 
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In the supergiant systems, the companion is of spectral type ear her than B2 and has evolved 
off the main sequence. Orbital periods are generally less than 10 days, orbits are circular, 
and the mass transfer takes place because of the strong stellar wind from the OB star 
and/or because of "incipient" Roche Lobe overflow. The Be star systems are characterized 
by emission lines (mainly the Balmer series) which originate in the equatorial circumstellar 
envelope of the companion star. The orbital periods in these systems tend to be longer than 
those of group a) and correlate well with the pulsar spin period (Corbet 1986). The orbits 
are usually moderately eccentric. Transient activity is common in Be HMXBs; different 
types of outburst have been observed from different sources, and occasionally also from the 
same source. Giant outbursts involve high peak luminosities, occur at any orbital phase and 
show only little (if any) orbital modulation of the X-ray flux. Recurrent outbursts usually 
involve lower peak luminosities, tend to occur close to periastron and sometimes show a 
strong modulation of the X-ray flux with the orbital phase (see e.g. Stella, White, & Rosncr 
1986). To date, about 80 HMXBs arc known; ~ 40 of them show periodic X-ray pulsations 
with spin periods distributed over a wide range from 69 ms to ~ 24 min (see Nagasc 1989 
for a review). 

While studying the atoll source GX 9-1-1 with BeppoSAX, we have discovered a new 
X-ray pulsator in a HMXB. We report here the results of the timing analysis of this new 
X-ray source. 

2. Observations 

The observation of the GX 9+1 fleld was carried out from 2001 September 16 02:01:30.0 
(UTC) to 2001 September 20 03:00:08.5 (UTC), using the co-aligned Narrow Field Instru- 
ments (NFIs) on board BeppoSAX. These are: a Low Energy Concentrator Spectrometer 
(LEGS; energy range 0.1-10 keV; Parmar et al. 1997), two Medium Energy Concentrator 
Spectrometers (MEGS; energy range 1-10 keV; Boella et al. 1997), a High Pressure Gas 
Scintillation Proportional Counter (HPGSPG; energy range 7-60 keV; Manzo et al. 1997), 
and a Phoswich Detector System (PDS; energy range 13-200 keV; Prontera et al. 1997). The 
exposure times were ~ 60 ks, ~ 149 ks, ~ 142 ks, ~ 71 ks, for LEGS, MEGS, HPGSPG and 
PDS, respectively. The circular fleld of view (FOV) of the LEGS and MEGS is 37 and 56 
in diameter, respectively, while those of the HPGSPG and PDS arc hexagonal with FWHM 
of 78' and 66', respectively. The LEGS and MEGS detectors arc position sensitive counters 
with imaging capability. The position reconstruction uncertainty for MEGS is 0.5' in the 
central area of 9' radius, and ~ 1.5' in the outer region of the FOV (Boella et al. 1997). The 
HPGSPG and PDS systems do not have imaging capabihties, and their data are therefore 
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difficult to interpret and analyse for individual sources when the FOV includes more than 
one source. 

Figure 1 shows the BeppoSAX/MECS image (1-10 keV) centered at the position of the 
bright atoll source GX 9+1. A fainter source is visible at R.A.(2000.0)= 18h 02m 39.9s and 
Dec. (2000.0)= —20° 17' 13.50" (position uncertainty 2'), at an angular distance from GX 
9+1 of ~ 22'. The source was outside the FOV of the BeppoSAX/LECS. We have verified 
its presence in both the MECS2 and MECS3 images, which probably excludes that this 
was a ghost image of a source outside the MEGS FOV. Moreover, we can be sure that the 
source was within the PDS FOV, because the source X-ray pulsations were also detected in 
the PDS data (see below). We searched for known X-ray sources in a circular region of 30' 
centered at GX 9+1 in the SIMBAD data base. We found no known sources with a position 
compatible with that of the faint source; we therefore designate this serendipitous source as 
SAX J1802.7-2017. 

We extracted the MEGS events of SAX J1802. 7-2017 from an eUiptical region centered at 
R.A.(2000.0)= ISh 02m 39.9s and Dec.(2000.0)= -20°17'l3.50" and the background events 
from an elliptical region similar to the one used for the source, centered at a symmetric 
position with respect to the center of the MEGS FOV and not contaminated by GX 9+1 
(see Fig. 1). In Figure 2 the lightcurves of SAX J1802. 7-2017 plus background (crosses) 
and of the background (stars) are shown. The source was not detected during the first 50 
ks when the count rate within the extraction region (0.034 + 0.001 c/s) was compatible with 
the background count rate. In the following ~ 300 ks the source count rate showed a large 
variability with an average intensity significantly above the background level. Two fiaring 
events took place ~ 110 ks and ~ 300 ks from the beginning of the observation. 



3. Temporal Analysis 

The arrival time of all events in the MEGS and PDS were corrected to the solar system 
barycenter. We searched for periodicities by computing a Power Spectrum Density (PSD) in 
the range between 4 x 10^'^ Hz and 1 kHz from Fast Fourier Transforms (FFTs) performed 
on MEGS data of SAX J1802.7-2017. 

A main peak at a frequency of ~ 7.20 mHz, corresponding to a period of ~ 139 s, is 
evident (see Fig. 3). The second harmonic is also clearly visible. We performed a folding 
search for periods centered around ~ 139 s finding a not well defined period, because two 
peaks were present at ~ 139.49 s and ~ 139.70 s, respectively. 

To study possible delays in the pulse arrival times, or, equivalently, pulse period vari- 
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Fig. 1.— BeppoSAX/MECS image (1-10 keV) of the GX 9+1 field. SAX J1802.7-2017 is 
visible at the top left. Its angular separation from GX 9+1 is ~ 22'. The extraction region 
of the source and background are also shown. The two semi-circular "cut-outs" are due to 
the removal of internal calibration source events. 




Fig. 2.— MEGS Lightcurves (1-10 keV) of SAX J1802. 7-2017 plus background (crosses) and 
of the background (stars). The bin time is 3000 s. 
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ations we divided the whole data set into 35 consecutive intervals each having a length of 
about 10 ks. We found that only 27 of these intervals had enough statistics to carry out the 
analysis described below. A folding search was performed in each interval for a range of trial 
periods centered around 139 s. The corresponding best periods were obtained by fitting the 
versus trial period curve with a Gaussian function. The best periods varied significantly 
and with continuity between 139.44 s and 139.86 s with an average value of 139.608 s. 

Phase delays (Fig. 4) were obtained by cross-correlating the folded lightcurves obtained 
for each of the 27 intervals with that of the whole observation used as a template (an average 
pulse period of 139.608 s was adopted in the folding). Initially we fitted the phase delays 

using 

A0 = ao + fll tn + 0.2 4, (1) 

where t„ is the arrival time of the nth pulse. In this formula the linear term ai = APp^/se/ Pp^i^^ 
and the quadratic term a2 = Ppuisef^Ppuise related to a correction to the pulse period and 
the derivative of the pulse period, respectively. We obtained a Xred 2.6 indicating that the 
modulation of the phase delays cannot only be explained by the presence of a derivative of 
the pulse period. Moreover we find a value of Ppuise/ Ppuise of 0.40 ± 0.01 yr~^, more than an 
order of magnitude larger than the largest measured value in any known X-ray pulsar (i.e. 
— 10~^ yr~^ in GX 1+4, see Pereira et al. 1999). We checked whether the modulation of the 
phase delays could be explained by the propagation delays due to the orbital motion of the 
X-ray pulsar around a companion star by fitting the phases with 



A(f) = Qq + tti tn + B COS 



'271 {tn -Tt,/2) 



(2) 



-'orb 

The results of the fit are reported in Table 1. The x^ed Since a.j, sini = Pp^iseB, we 

find sini ~ 70 It-s; the corrected pulse period was ~ 139.612 s. 

We tried to add a quadratic term to expression (2) to take into account a spin period 
derivative, but, because of the small number of points and a relatively short observation, we 
could not constrain all parameters of the fit. In any case, by fixing the quadratic term in the 
range —9 x 10~^ - 9 x 10~^ days~^ (corresponding to a Ppuise/ Ppuise in the range — 10~^ - 
10~^ yr~^) we found that the orbital parameters do not change significantly (~ 2% changes 
in the orbital period value and ~ 4% in the ax sini value) with respect to the parameters 
reported in Table 1. 

We corrected the arrival times of all the events, observed by the MEGS and PDS in- 
struments, to the center of mass of the binary system using the orbital parameters and the 
corrected pulse period reported in Table 1. We then re-computed, on the corrected folded 
lightcurves, the phase delays, and found that they are compatible with zero phase delay 
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Fig. 3.— Power spectrum density of SAX J1802.7-2017 obtained from MECS data (1.0-10.5 
keV energy band) extracted from an elliptical region centered on the source as described in 
the text. 61 PSDs, computed from time intervals of ~ 5 ks length, were averaged. The 
fundamental and second harmonic are clearly visible at 7.2 mHz and 14.4 mHz, respectively. 
The PSD is normalized according to prescription of Leahy et al. (1983). 
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Fig. 4. — Phase delays of SAX J1802.7-2017 as a function of time. The sohd hne is the best 
fit function (2) discussed in the text; the orbital period is ~ 4.6 days. 
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Table 1: Orbital parameters of SAX J1802. 7-2017. Errors are at 1 cr confidence level. 



Parameter Value 

Vo 0:0431 

ai 0.02 ± 0.03 days"^ 

B 0.501°:°^ 

Porb 4.61°:^ days 

Ux sin i 70^7° It-s 

52168.22l°i° MJD 

P 1 "^Q R1 9+0.006 

^pulse -'-•^"^■"-'-^-0.007 ° 



(see Fig. 5). This suggests that our data are compatible with a circular orbit, although an 
eccentricity could be present but not appreciated because of the relatively low statistics of 
our data. To estimate a rough upper limit to the eccentricity of the system, we fitted our 
phase delays substituting the circular orbital correction in equation (2) with a first-order 
approximation of the eccentric orbital correction (see e.g. van der Klis & Bonnet-Bidaud 
1984), thus obtaining e < 0.2 (90% confidence level) and xfcd ~ 

Using the estimated pulse period we folded the MEGS lightcurves in the energy bands 
1-3 keV, 3-6 keV and 6-10 kcV, and the PDS lightcurves in the energy bands 13-25 keV and 
25-80 keV, excluding the first 50 ks of our observation. The main pulse of the modulation 
presents a peak around phase 0.6. A secondary peak is visible around phase 0.1, which 
appears to be more prominent at higher energies. The pulse fraction, defined as {Imax ~ 
Imin)/ I-max-i with I^ax ^ud Imin the maximum and minimum count rate, are (43 ± 7)%, 
(57 ± 6)%, (62 ± 8)%, (7 ± 2)%, and (36 ± 3)% in the 1-3 keV, 3-6 keV, 6-10 keV, 13-25 
keV, and 25-80 keV energy band, respectively. The pulse fraction in the 13-25 kcV and 
25-80 keV energy bands is strongly reduced by the presence of GX 9+1 in the PDS FOV 
(note however that the pulse fraction between 25-80 keV is somewhat higher because GX 
9+1 is weaker at these energies). In Figure 6 we show the folded lightcurves in the energy 
bands 1-3 keV, 3-6 keV, 6-10 keV and 25-80 keV. 



4. Discussion 

We have discovered a serendipitous source, SAX J1802. 7-2017, in a BeppoSAX/NFI 
observation of the bright atoll source GX 9+1. SAX J1802. 7-2017 shows coherent pulsations 
at a period of 139.612 s, indicating that the compact object in this system is most likely 
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Fig. 5. — Phase delays after the barycentric correction with respect to the center of mass of 
the binary system. 
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Fig. 6. — Folded lightcurves in four different energy bands. From the top to the bottom: 1-3 
keV, 3-6 keV, 6-10 keV (MEGS), 25-80 keV (PDS). The folding is obtained for a pulse period 
of 139.612 s and the zero epoch was assumed at the superior conjunction Tt^i2 — 52168.22 
MJD. 



9 



an accreting magnetic NS. Pulse arrival times show a sinusoidal modulation at a period of 
4.6 days, which wc interpret as due to the orbital motion of the source. The values of the 
orbital period and of the pulse period, as well as the pulse fraction in the MECS energy 
range (~ 40% — 60%), arc consistent with the values usually found for HMXB systems (e.g. 
Oosterbroek et al. 1999; Bildsten et al. 1997). 

The mass function of the system is 



where Mc is the companion star mass, G is the gravitational constant and q is the ratio 
of the compact object mass to the companion mass. The value of the mass function is not 
compatible with the system being a cataclysmic variable in which a white dwarf is orbiting 
a ~ 1 Mq companion. For a NS mass of 1.4 M© we have estimated a lower hmit on the 
companion star mass, Mc > 11 Mq (90% confidence level), typical of a HMXB. The epoch 
of the NS superior conjunction (which would correspond to the mid-eclipse time) derived 
from the fit of the orbital parameters falls during the first 50 ks of our observation, when 
the count rate of SAX J1802. 7-2017 is compatible with the background count rate. This is 
compatible with the possible presence of an eclipse during the first 50 ks. Assuming that 
the eclipse ends ~ 0.58 days after the beginning of our observation (52168.69 MJD, sec Fig. 
2) and taking into account that the zero epoch is 0.11 ± 0.10 days, the half duration of the 
eclipse should be 0.47 ± 0.10 days. The half-angle of the echpse subtending the portion of 
circular orbit covered by the source is 9 — 0.64 ±0.14 rad. The eclipse duration is related 
to the inclination angle i of the system by: Rc/a — (cos^i + sin^ i sin^ ^)^/^, where Rc is 
the radius of a spherical primary star and a is (approximately) the separation between the 
two stars (see Primini et al. 1976). Then a rough lower limit on the companion star radius 



To estimate the luminosity of the source we performed a spectral analysis in the MECS 
range using the correct effective area for the position of the soTirce in the MECS FOV (G. 
Cusumano, private communication). Although the low statistics of the data did not allow an 
accurate spectral analysis, we inferred a luminosity of ~ 5.6 x 10^^ erg/s over the 1.8-10 keV 
range by fitting the data to a power law (photon index —0.1 ±0.1 and Nh = x 10^^ 

cm~^) and assuming a source distance of 10 kpc. The lightcurve of the system shows a 
large variabihty suggesting that the system could be accreting via a stellar wind. This is 
also suggested by its long pulse period and relatively short orbital period, typical of HMXBs 
accreting via stellar wind such as Vela X-1 and 4U 1538-52 (see e.g. Corbet 1986; Nagase 
1989). Indeed SAX J1802.7-2017 shows similarities to 4U 1538-52, which has a long spin 
period of ~ 529 s, a short orbital period of ~ 3.7 days, and a mass of 19.8 ± 3.3 Mq 




is Rc > 14 Rq. 
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for the companion star (Reynolds, Bell, & Hilditch 1992). However, while the averaged 
luminosity of 4U 1538-52 and Vela X-1 is around 5 x 10^^ ergs s~^, the averaged luminosity 
of SAX J1802.7-2017 measured during our observation is at least one order of magnitude 
lower (assuming an upper limit to the distance of 10 kpc). The source seems therefore to 
be quite underluminous with respect to massive binaries accreting via stellar wind, although 
more observations are needed to confirm this conclusion. 

Note that, as the source has not been reported by other X-ray astronomy satellites 
before (in particular we did not find any detected source at the position of SAX J1802.7- 
2017 in ROSAT images, although the field of GX 9+1 was not observed by ASCA) it might 
well be a transient HMXB. Most of the transient HMXBs are Be-star systems (Liu et al. 
2000). However, if SAX J1802. 7-2017 belongs to this class it should be atypical considering 
that its orbital period and pulse period do not follow the Corbet correlation (Corbet 1986). 
Note also that the presence of an eccentricity cannot be excluded by our data, and more 
observations are needed to confirm the orbital parameters of the source. 

We thank Dr. Giancarlo Cusumano for providing us the corrected effective area for our 
source and the anonymous referee for useful suggestions. This work was partially supported 
by the Italian Space Agency (ASI), by the Ministero della Istruzione, della Universita e della 
Ricerca (MIUR) and the Netherlands Organization for Scientific Research (NWO). 
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